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PACS 32 . 80 . Lg - Mechanical effects of light on atoms, molecules, and ions 

PACS 05 . 40 . Jc - Brownian motion 

PACS 32 . 80 . Pj - Optical cooling of atoms; trapping 

Abstract. - We study the influence of the lattice topography and the coupling between motion 
in different directions, for a three-dimensional Brownian motor based on cold atoms in a double 
optical lattice. Due to controllable relative spatial phases between the lattices, our Brownian 
motor can induce drifts in arbitrary directions. Since the lattices couple the different directions, 
the relation between the phase shifts and the directionality of the induced drift is non trivial. Here 
is therefore this relation investigated experimentally by systematically varying the relative spatial 
phase in two dimensions, while monitoring the vertically induced drift and the temperature. A 
relative spatial phase range of 27r x 27r is covered. We show that a drift, controllable both in 
speed and direction, can bo achieved, by varying the phase both parallel and perpendicular to the 
direction of the measured induced drift. The experimental results are qualitatively reproduced by 
numerical simulations of a simplified, classical model of the system. 



tZ3 
>> 

Introduction. — A Brownian motor (BM) converts 
I ^f andom fluctuations into deterministic work [1-3]. Brow- 
nian motors exits naturally as, e.g., protein motors and 
CN intra-cell motion [4] , and a general understanding of their 
^mechanism is of fundamental interest. The rectification 
^ mechanism requires that the system is both brought out of 
thermal equilibrium [1] and spatially or temporally asym- 
metric [1,5]. Although it has not been proven, fulfilling 
these two requirements arc generally sufficient to realise a 
BM. These requirements can be met using ultra cold atoms 
trapped in optical lattices [6,7]. Generally, the symmetry 
is broken by using a spatially asymmetric or flashing po- 
^ tential. In this type of noise rectifier, the direction of the 
induced drift is often fixed for a given potential or control- 
KNlablc in just ID [1-3,8-12]. 

^ Our BM [13, 14] is based on two symmetric potentials, 
with an asymmetry that originates from a combination of 
a relative spatial phase shift between the potentials and an 
unequal transfer rate between them. Our system possesses 
an inherent ability to induce drifts in an arbitrary direc- 
tion in three dimensions for the ultra cold atoms interact- 
ing with the potentials. The direction is mainly controlled 
by the relative spatial phase between the two potentials. 
In [13,14] we demonstrated that our BM does work in more 
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than one dimension, but how this works was not investi- 
gated. This question is especially complicated since the 
four-beam configuration induces a coupling between the 
different spatial dimensions. This coupling affects the 3D 
behaviour of our BM profoundly and is the key to any un- 
derstanding and control of the three-dimensional aspects 
of our BM. 

In this paper, we present a study of the dimensional 
coupling between the relative spatial phases and its in- 
fluence on the induced drift, in two dimensions. This is 
done both experimentally and with numerical simulations. 
This study results in a better understanding of the mul- 
tidimcnsionality of our BM, and in a good control to set 
the induced drift to a chosen speed and direction. This 
also renders a controlled dynamical BM possible, where 
a time-dependent phase can induce realtime drifts along 
any pre-chosen trajectory. It even allows to control the in- 
duced drift in the vertical direction by purely controlling 
the phase shift horizontally, and vice versa. By studying 
the phase dependence, we can also study the coupling be- 
tween the dimensions of the two potentials. Information 
about what potentials the atoms really experience and the 
adiabaticy [7] of the interaction can therefore be investi- 
gated. The numerical simulations are done with a simpli- 
fied, classical model, which mimics the main features of 
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our BM. This qualitatively reproduce the main features of 
the results from the experiments. 

System. — The potentials used are realised from a 
double optical lattice (DOL) [15, 16], interacting with 
^^■^Cs. An optical lattice is a periodic potential where 
atoms can be cooled and trapped due to dissipation and 
light shifts [6], and is created from the interference be- 
tween two or more light fields. A DOL consists of two spa- 
tially overlapped state-dependent optical lattices [15, 16]. 




Fig. 1: (Colour online) Simplified, classical model of the BM 
mechanism, (a) (/p = 0, the system is spatially symmetric and 
the particle undergoes Hamiltonian oscillation at the bottom of 
the wells. No drift is induced, (b) if ^ 0, the spatial-temporal 
symmetry is broken and a drift v to the left is induced. In both 
cases there is a large difference in the transfer rates 7 between 
the lattices, indicated with arrow of different thickness in the 
figure. Note that this a strongly simplified model, presented to 
provide an understanding of the phase dependence. Features 
such as diffusion and friction, vital for the BM, are not present 
in this model. 

To qualitatively understand the induced drift depen- 
dence on the relative spatial phase, a simple classical 
model is used [17, 18], see figure 1. Consider a classi- 
cal Brownian particle situated in either of two symmetric 
and periodic potentials (C/a and Ub), coupled with un- 
equal transfer rates (7a^b and 7b^a)- The particle will 
shift back and forth between the potentials, but on average 
spend longer time in one of them. When the potentials are 
in phase, the trapped particle will jump between potentials 
and undergo Hamiltonian oscillations near the bottom of 
them, see figure la. If the potentials are given a non-zero 
relative spatial phase shift, see figure lb, with transfer 
rates properly chosen, the particle will spend long enough 
time in the long-lived potential in order to be located close 
to its minimum. After a certain time the particle is trans- 
ferred to the short-lived potential, where it will be located 
at a slope of the potential. After a short while, it will 
be pumped back to the long-lived potential. The particle 
will therefore, on average gain kinetic energy in a certain 
direction when jumping between the potentials [18]. The 



optimal ratio between the transfer rates has been exper- 
imentally determined to 7a^b : 7b^a — 9:1 [14]. This 
is also related to the time atoms spend in each potential. 
The time spent in potential A should optimally be of the 
order of the the inverse of this potential's oscillation fre- 
quency [19], while the optimal time spent in potential B 
is much larger. The oscillation frequencies are typically of 
the order of 100 kHz. 

Methodology. — Detailed descriptions of the general 
experimental setup are given in [15,16,20,21]. More details 
concerning the construction and the control of the BM are 
found in [13, 14,22]. In short, cesium atoms are first accu- 
mulated in a magneto-optical trap, from where they are 
transfered to the double optical lattice. Both lattices (A 
and B) have a 3D topography and are constructed from 
four beam configurations, in which two beams are prop- 
agating in the a;2-plane, and two in yz-plane. All beams 
have a 45 ° angle with respect to the vertical 2-axis with 
polarizations perpendicular to the plane of propagation 
("3D lin_Llin configuration" [6]). The lattices are state- 
dependent and operate from different hyperfine ground 
states within the Cs D2 line (6s ^Si/2 6p ^P3/2). Lat- 
tice A operates close to the i^g = 3 Fo = 4 transition 
and lattice B operates close to the Fg — 4 = 5 

transition. The irradiances and the frequencies of the four 
beams are chosen to generate pumping rates between the 
lattices that maximise the induced drift. 

The optical path lengths are controllable in all beams, 
resulting in a control of the relative spatial phase. This de- 
pendence originates from a slight difference in wavelength 
between the lattices. The difference is small enough for the 
relative spatial phase to be effectively the same through- 
out the whole atomic cloud [15,22]. 

The induced drift velocity is obtained by a ballistic time- 
of- flight method (TOF), where the trapped atoms are re- 
leased by quickly turning off the lattice beams. The atoms 
then freely fall through a resonant probe located approxi- 
mately 5 cm below the trapping region. The fluorescence 
from the atoms passing through the probe is detected and 
gives a measurement of the arrival time t of the atomic 
cloud [23,24]. Since the drift velocity is constant [13], the 
arrival time can be converted into drift velocity through 



where g is the gravitational constant, r is the interaction 
time of the atoms with the DOL and Sp, is the distance 
to the probe. 

Results. — To investigate the relative spatial phase 
dependency experimentally, we change the phase along z 
{^Pz) and y {(py) by slightly more than 27r. This is done 
in about 20 steps in each direction, covering all steps in z 
for each y value. Here, ipz^y = corresponds to a perfect 
overlap of the two lattices. Since these intervals exceed 
2tt, whole periods along both the z and y axes are covered. 
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For each combination phases, the TOF signal is measured. 
To improve statistics, an average over five TOF measure- 
ments is made for each combination. A compilation of 
the measurements is shown in figure 2. Also displayed is 
the temperature measured at each step, which is used to 
determine the phase. 
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Fig. 2: (Colour online), (a) Measured induced drift, and (b) 
temperature, as a function of the relative spatial phase along 
the y- and 2;-axes. A difference in the pattern created by the 
induced drift maxima and minima is evident. The dashed lines 
indicates the cuts shown in figure 3. 

The origin, where (py^z is zero, is chosen at a minimum 
of the temperature, since a perfect spatial overlap of the 
lattices is assumed to correspond to a temperature min- 
imum [15, 16]. The scale is obtained by measuring the 
periodicity of the temperature and comparing it to the 
change of the optical path lengths. The measurement in- 
dicates that the temperature maxima occur when the di- 
abatic potentials [6, 7] are shifted with a relative spatial 
phase of tt. Since the induced drift velocities and the tem- 
perature are measured simultaneously, the relative phase 
scale obtained from the temperature is applied to the BM. 

Throughout the measurement, lattice B is red-detuned 
below resonance by 40r and has an irradiance of 0.2 



mW/cm^ in each beam and lattice A is red-detuned by SSF 
and has an irradiance of 0.7 mW/cm^ in each beam, where 
F is the natural line width of the excited state (F/27r — 5.2 
MHz). These parameter values are chosen to optimize the 
induced drift [14] while still having a strong TOF signal 
with a good signal to noise ratio. 
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Fig. 3: Cuts in the two dimensional plot of figure 2a. (a) Cut 
along j/-axis, at ~ ^tt. A clear splitting of the minima is 
visible at ^Sj^ ~ tt (b) Cut along z-axis, at ~ 0. The induced 
drift shows a close to sinusoidal shape, (c) A diagonal cut, at 
(pz ~ A doubled frequency is evident. 
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Even though the induced drift has the same periodicity, 
the topography along the z and y axes differ. Along the z 
axis, the drift has a close to sinusoidal dependence on ipz, 
for any choice of (py (see figure 3b) , although the amplitude 
and the phase of the drift velocity strongly depend on the 
relative spatial phase along the y axis. 

Along the y axis, however, the induced drift clearly de- 
viates from a sinusoidal shape, having two minima per 
period (see fig 3a) . Such a difference between the z and y 
intersections is expected since the lattice structure differs 
slightly along these directions [6]. By moving diagonally 
in the yz phase plane, the periodicity is doubled (see fig 
3(c)). 

In order to understand the qualitative behaviour of our 
Brownian motor, we have performed simulations, using a 
simplified, classical model of our system [17, 18,25]. The 
two-dimensional numerical simulations are done for a clas- 
sical Brownian particle situated in either of two potentials, 
corresponding to the lowest diabatic potential [15, 17] of 
each lattice in a DOL configuration, identical to the one 
used in the experiments. The dynamics of the atom is 
obtained from the Fokker-Planck equation [26,27]. The 
induced drifts in the z and y directions as a function of 
the relative spatial phase are shown in figure 4. The simu- 
lations show a clear periodic behaviour in both the z and 
y directions, and the fact that the induced drift maxima 
in positive and negative directions do not lie on a straight 
line along the y axis shows that the potentials in differ- 
ent directions are coupled. That is, it is not possible to 
go from a maxima to a minima by purely changing (py, 
ipz must vary as well. This is also visible in the equation 
describing the potentials, 

=F 44cos(fcj;2:) cos{kyy) cos(/czz)}, (2) 

where Aq is the scattering rate and k is the effective an- 
gular wave vectors along the axes [6]. This equation also 
shows that the potentials are identical in x and y, which 
would give identical relative spatial phase dependencies 
in X and y. The simulations agree with the measurement, 
where a clear spatial coupling between the potentials is ev- 
ident. Both also show a clear pattern, and even though the 
pattern does differ, this confirms that we can control our 
BM via the phase in arbitrary directions. While an exten- 
sion of the simulations to 3D is conceptually straightfor- 
ward, it would be computationally demanding. The third 
dimension would, if anything, just enhance the coupling 
observed in 2D. Moreover, while this model contains the 
basic ingredients of our BM (diffusion, friction, coupled 
potentials, etc.), it neglects potentially important factors 
(such as the spatial dependence of diffusion and friction or 
the presence of a manifold of potentials due to the mag- 
netic substates), so no far reaching conclusions should be 
drawn from the comparison of the simulations with the 
experimental results. 
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Fig. 4: (Colour online) Simulation of the induced drift as a 
function of the relative spatial phases ifz and ipy. The speed is 
given in mm/s. (a) False colour contour plot of induced drift in 
z. (b) Arrow plot of the induced drift in z and y. The arrows 
indicate the size and the direction of the induced drift. The 
phase at wich the velocity reaches its maximum differs along 
the z and y directions. This agrees with the result obtained 
in [14]. 



Discussion. — The induced drift of our BM can be 
controlled by the difference between the transfer rates 
between the lattices, the diffusion, the potential height 
and the relative spatial phase between the lattices. The 
dependence of the first three parameters is investigated 
in [13, 14, 25]. The phase dependency was there fully in- 
vestigated for drifts measured along only in one dimension, 
and phases varied along the same direction. Drifts in more 
than one dimension were demonstrated but drifts induced 
by a phase shift in a perpendicular direction were not in- 
vestigated at all. Since the potentials couples the different 
directions, the non-trivial relation between the phase shift 
of the lattices and the directionality of the induced drift is 
at the heart of the multidimensionality of our BM. Such 
an investigation is therefore of fundamental interest and 
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necessary obtain any understanding or control of our BM. 

By surveying the vertically induced drift by 27r x 2-k 
in the yz phase plane, the influence of the lattice topog- 
raphy on the structure and periodicity of the drift can 
be determined. Due to the identical potentials in the x 
and y directions, a generalisation to three dimensions is 
possible. The induced drift dependence on the transfer 
rates, i.e., the irradiances and detunings, has earlier been 
investigated [13,14], and thus we can now fully describe 
our three-dimensional BM, and with good control set the 
induced drift velocity to any arbitrary direction, with a 
controllable magnitude in three dimensions. Such control 
is necessary in future experiments where fast phase shifts 
should render it possible to realise a dynamical Brownian 
motor, where the atoms can change direction in real time 
and therefore travel along more or less any pre-chosen tra- 
jectory, limited only by the lifetime of the DOL-caused dif- 
fusion. These fast phase shifts can be realised with electro- 
optical modulators (EOM). The implementation of these 
EOMs also makes it possible to scan the phase over a range 
of 2ti in all three dimensions simultaneously, which is not 
possible with the current set-up. 

The multidimensional coupling between the phase de- 
pendency of the lattice topography and the induced drift 
render it possible to redirect the motion in a certain di- 
rection without changing the phase in that direction. The 
induced drift can, by changing the phase in an orthogonal 
direction, be enhanced, put to zero or even be inverted. 
This is clearly evident in flgure 2(a), and has not been 
shown before. 

The measured periodicity of the temperature also con- 
flrms that the diabatic potentials are the relevant poten- 
tials in the system, as shown in [14]. The two potentials 
in the classical simulations are therefore chosen as the two 
lowest diabatic potentials of the two state-dependent po- 
tentials of the DOL used in the experiment [17]. The sim- 
ulation qualitatively reproduces the main features of our 
BM although it is based on a classical model. A closer 
comparison with the experimental result show clear dif- 
ferences, especially the location of the maxima in induced 
drift velocity in the z direction. This shows that the simple 
classical picture is inadequate for extracting more detailed 
information about our system, and that not just the lowest 
diabatic potential is relevant. In reality an optical lattice 
works on a manifold of potentials, each corresponding to a 
transition between different magnetic sub-states [7]. This 
manifold of potentials is also of importance in the model 
explaining the relative spatial phase dependence of the 
temperature in a DOL [16], seen in flgure 2b. The depen- 
dence of the temperature on the phase also affects our BM, 
since a change in temperature reflects different diffusion 
and friction in the system. This coupling may indeed con- 
tribute to the double peak structure of the maxima in the 
induced drift in the negative z-direction, since the split- 
ting overlaps with the temperature maxima. Moreover, 
the phase in the x direction is roughly set to zero at the 
beginning of the experiment by measuring the minimum 



in temperature. During the scan this value may change 
slightly but the effect is small, as long as the phase is small. 
Finally, a full physical understanding of the structure of 
our BM calls for a quantum model where the entire man- 
ifold of potentials that exist in the real DOL is accounted 
for. Such a model is under construction. 

Conclusion. — In summary, we have shown that the 
relative spatial phase of the lattices, due the the four beam 
conflguration of the latter, couple dimensionally the in- 
duced drift. This affects the multidimensional behaviour 
of our Brownian motor profoundly. These effects have 
been investigated and we showed that the vertical induced 
drift can be controlled by either of the vertical or the hor- 
izontal phase, or a combination of the two. This was done 
by measuring the vertical drift velocity in a 27r x 27r area in 
the yz phase plane. Prom the measurement, the structure 
and periodicity of the drift was extracted. This gives us 
full control over our three dimensional Brownian motor. 
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